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Abstract 

A new point of view about the deep origin of thermodynamic phase 
transitions is sketched. The main idea is to link the appearance of phase 
transitions to some major topology change of suitable submanifolds of 
phase space instead of linking them to non-analyticity, as is usual in 
the Yang-Lee and in the Dobrushin-Ruelle-Sinai theories. In the new 
framework a new possibility appears to properly define a mathematical 
counterpart of phase transitions also at finite number of degrees of free- 
dom. This is of prospective interest to the study of those systems that 
challenge the conventional approaches, as is the case of phase transitions 
in nuclear clusters. 

1 Introduction 

The current mathematical definitions of thermodynamical phase tran- 
sitions are based on the loss of analyticity of thermodynamical observ- 
ables. This is to some extent suggested, though not implied, by the 
experimental relations among macroscopic observables. This conflicts 
with the analytic character of the statistical weights that are attributed 
by any ensemble in statistical mechanics to the microscopic configura- 
tions. Thus, as it has been proved by the Onsager solution of the 2D 
Ising model and by the Yang-Lee theorem, the only way to elimiwate 
this conflict is to consider the limit N oo (thermodynamic limit) □. 

Obviously phase transitions in Nature occur at finite A'^, but it is 
commonly argued that - for macroscopic objects - is so large that 
it can be considered "infinite" from a physical point of view. However, 
by looking at a small and embroidery-like snowflake that melts into a 
drop of water one can wonder why such a phenomenon should ever be 
explained only in terms of the infinite N limit. 

Moreover, there is a growing experimental evidence that phase tran- 
sitions can occur also at finite and small N, i.e. with N -C Avogadro 
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number. Some examples of small TV systems undergoing phase transi- 
tions are: i) nuclear clusters as well as atomic and molecular clusters; 
ii) nano and mesoscopic systems; in) polymers and proteins; iv) small 
drops of quantum fluids (BEC, superfluids and superconductors). 

Whence the prospective interest of a new mathematical characteri- 
zation of thermodynamical phase transitions which, instead of resorting 
to the loss of analyticity of macroscopic observables, might naturally 
encompass also finite A'' systems. The search for a broader mathemat- 
ical definition of phase transitions is also of potential interest to the 
treatment of other important topics in statistical physics, as is the case 
of amorphous and disordered systems (like glasses and spin-glasses) , or 
for a better understanding of phase transitions in the microcanonical 
ensemble, for first-order phase transitions, and so on. 

1.1 Heuristic arguments 

Everything here refers to classical Hamiltonian systems with continuous 
variables and described by a standard Hamiltonian 

H[p= {pi,. . .,pN),q = {qi,-. ■ ,qN)] = '^^pl + V{q) . (1) 

i 

As a consequence of a systematic study of the dynamical counterpart of 
thermodynamical phase transitions, performed by numerically solving 
the Ifamilton equations of motion of (jl|), it has been found that Lya- 
punov exponents dispJay "singular" energy and temperature patterns at 
the transition point 0. Thus, since in a differential-geometric descrip- 
tion of Hamiltonian chaos Lyapunov exponents are seen as probes of 
the geometry of certain submanifolds of configuration space, it has been 
conjectured that a phase transition could be due to some major geomet- 
rical, ancUBOKsibly topological, change in the support of the statistical 
measures D'DU. This is to say that it is on the basis of the just mentioned 
work that we have been led to formulate the following argument. 
Let us consider the canonical configurational partition function 

where a co-area formula has been used to unfold the structure integrals 

i.e an infinite collection of purely geometric integrals on E^~^, the 
eqmpotential hypersurfaces of configuration space defined by 'S^'^ = 
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If we consider the microcanonical ensemble, the basic mathematical 
object is the phase space volume 



where 



n^-\E-^)= I d''qe[V{q)^{E^rj)]= I dvl ^(4) 



IVVI 



whence 



(2^ f _Ja_ 

vnf) Jo L iivvii 



n{E) = / dv '-^^^ I dv I ^ , (5) 



also here, as in the above decomposition of Zc{l3, N), the only non-trivial 
objects are the structure integrals (^. 

Once the microscopic interaction potential V{q) is given, the con- 
figuration space of the system is automatically foliated into the family 
{E„}„g27j of equipotential hypersurfaces independently of any statistical 
measure we may wish to use. Now, from standard statistical mechani- 
cal arguments we know that the larger is the number of particles the 
closer to some E„ are the microstates that significantly contribute to the 
statistical averages of thermodynamic observables. At large A'', and at 
any given value of the inverse temperature /3, the effective support of 
the canonical measure is narrowed very close to a single E„ — E„{^^); 
similarly, in the microcanonical ensemble, the fluctuations of potential 
and kinetic energies tend to vanish at increasing A'' so that the effective 
contributions to i^{E) come from a close neighborhood of a E„ = E„(£;^)- 

Now, the "topological conjecture" consists in assuming that some suit- 
able change of the topology of the {E„}, occurring at some Vc = Wc(/3c) 
(or Vc = Vc{Ec)), is the deep origin of the singular behavior of thermody- 
namic observables at a phase transition; (by change of topology we mean 
that {E„}„<„^ are not diffeomorphic to the {E„}„>„^). In other words, 
the claim is that the canonical and microcanonical measures must "feel" 
a big and sudden change - if any - of the topology of the equipotential 
hypersurfaces of their underlying supports, the consequence being the 
appearance of the typical signals of a phase transition, i.e. almost sin- 
gular energy or temperature dependences of the averages of appropriate 
observables. The larger is A'' the narrower is the effective support of the 
measure and hence the sharper can be the mentioned signals. Eventually, 
in the N —> oo limit this sharpening will lead to non-analyticity. 
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2 A theorem about topology and phase transitions 

We have recently proved a theorem stating that topological changes of 
hypersurfaces E„ are necessarily at the origin of phase transitions 
It applies to physical systems described by short-range potentials V, 
bounded below, of the general form: 

V{{q}) = Yl yoiUc ^ q,\\) • 

(a, 7) a 

The theorem is proved in the reversed formulation: if the surfaces E,, 
with V = V/N £ I — \v^,Vm] are diffeomorphic, then no phase transi- 
iion will occur in the corresponding temperature interval [l3{v,n), P{^m)]- 
The proof is lenghty and rather complicated but it proceeds along a log- 
ically simple path. Diffeomorphicity of the Ei,, after the "non-critical 
neck theorem" in Morse theory, implies the absence of critical points of 
V, i.e. VV 7^ 0. In the absence of Morse critical points: 



dv'^ [J^^ \\vv\\j j^^-^ \ \\vv\ 



da 



where 



-D^(IIVVir') = 2\\VV\\~'^Mi - \\VV\\~^AV 

and where Mi is the trace of the shape operator of the hypersurface Et, . 
There is an operator algebra to generate the powers D'^ so that, being 
SN{if) ~ j^IuQn iv{v-)) the microcanonical configurational entropy, it 
is possible to show that 

sup SNi^)<oo ; sup ^^^^{v) < 00 ,k — 1, . . . ,4 

N,vei N,vei 

that is: if the E„ are diffeomorphic then Sn{^) is uniformely convergent 
in C^{I) as N ^ 00, and, by using the Legendre transform relationship 
<5']v('») = fN{/3) + P ■ V- + o{N) between microcanonical configurational 
entropy and the canonical free energy fniP) = -^\n Zc{(3, N), this im- 
plies that - in the N —> cx> limit - the canonical configurational free 
energy is /oo(/3) £ C^(/), i.e. at least twice differentiable, thus there 
are neither first nor second order phase transitions according to their 
standard definition. There is not a one-to-one correspondence between 
phase transitions and topology changes of the E„, the latters are neces- 
sary but not sufficient. Sufficiency conditions, to point out the special 
class of topolog¥,|Ghanges that give rise to phase transitions, are based 
on relations likeuO 



dnjv(u) f Ml da ,/ 0(^-1) \]1/n 
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that bridge thermodynamics and topology; here are the Betti 

numbers (cohomological invariants) of S„. It turns out that a '^second 
ordej" topology change, i.e. a sudden change in the way of changing of 
topology as a function of v, is sufficient to entail a first or a second order 
phase transition. Such topological discontinuities can exist and can be 
detected at any finite N, though they yield non-analyticity of thermo- 
dynamic observables only when the support of the statistical measure 
indefinitely shrinks with A*' oo. In other words, we have here the 
possibility of properly defining phase transitions also at finite A*', whose 
detection, direct or indirect, can be performed through quantities that 
probe the topology of the E„ . It is an interesting and open question how 
to make the link with other approaches tackling the finite riV transitions 
in a macroscopic phase space of thermodynamic variables u. 

3 A direct numerical confirmation 

The family of {'Sv}veiiz associated with ft V'* model on a d-dimensional 
lattices Z'^ with d — 1,2 has been usedu for a numerical check of the 
scenario sketched in the previous section. The model is described by 

(-T'^? + i'^')+ E (6) 

(ik) stands for nearest-neighbor sites, and in d = 2 it undergoes a second 
order phase transition. In order to directly probe if and how the topology 
change - in the sense of a breaking of diffeomorphicity of the surfaces 
E„ - is actually the counterpart of a phase transition, a diffeomorphism 
invariant has to be computed. This is a very challenging task because 
of the high dimensionality of the manifolds involved. One possibility 
is afforded by the Gauss-Bonnet-Hopf theorem that relates the Euler 
characteristic x(^u) with the total Gauss-Kronecker curvature Kg of 
the manifold 

X(E„) =7 j Kg da (7) 

which is valid in general for even dimensional hypersurfaces of euclidean 
spaces TZ^ , and where 7 = 2/l^oZ(S") is twice the inverse of the volume 
of an n-dimensional sphere of unit radius, and da is the invariant volume 
measure induced from IJZ^ . In Fig. ^ x(^") is reported vs v: the Id case 
gives a "smooth" pattern of whereas the 2d case yields a cusp- 

like shaped x('^) the phase transition point. There is here a direct 
evidence of a major and very sharp '^second order" topological transition 
that underlies the phase transition, it is also remarkable that with a 
very small lattice of A'' = 7 x 7 sites such a sharp signal-jwould never be 
obtained through standard thermodynamic observablesu. 
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Figure 1; Euler characteristic xi^v) for 1-d and 2-d ip^ lattice models. Open circles: 1-d 
case, N = 49; full circles: 2-d case, N = 7 x7. The vertical dotted line, computed separately 
for larger N, accurately locates the phase transition. Data are from Ref . [8] . 



References 

1. C.N. Yang and T.D. Lcc, Phys. Rev. 87, 404 (1952); D. Ruelle, 
Thermodynamic formalism, Encyclopaedia of Mathematics and its 
Applications, (Addison- Wesley, New York, 1978). 

2. L. Caiani, L. Casetti, C. Clementi, and M. Pettini, Phys. Rev. Lett. 
79, 4361 (1997); L. Caiani, L. Casetti, C. Clementi, G. Pettini, M. 
Pettini, and R. Gatto, Phys. Rev. E 57, 3886 (1998); M. Corruti- 
Sola, C. Clementi and M. Pettini, Phys. Rev. E 61, 5171 (2000). 

3. R. Pranzosi, L. Casetti, L.Spinelli and M. Pettini, Phys. Rev. E 
60, R5009 (1999); L. Casetti, E.G.D. Cohen and M. Pettini, Phys. 
Rev. Lett. 82, 4160 (1999). 

4. L. Casetti, E.G.D. Cohen and M. Pettini, Phys. Rep. 337, 237-342 
(2000). 

5. R. Pranzosi, (PhD thesis, Universita di Firenze, 1998). 

6. L. Spinelli, (PhD thesis, Universite de Provence, 1999). 

7. D. H. E. Gross, Phys. Rep. 279, 119 (1997). 

8. R. Pranzosi, M. Pettini and L. Spinelli, Phys. Rev. Lett. 84, 2774 
(2000). 

9. L. Caiani, L. Casetti, and M. Pettini, J. Phys. A: Math. Gen. 31, 
3357 (1998). 



6 



